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Abstract 

We investigate a circular cavity billiard within which a pair of identical hard disks of smaller 
but finite size is confined. Each disk shows a free motion except when bouncing elastically with its 
partner and with the boundary wall. Despite its circular symmetry, this system is nonintegrable and 
almost chaotic because of the (short-range) interaction between the disks. We quantize the system 
by incorporating the excluded volume effect for the wavefunction. Eigenvalues and eigenfunctions 
are obtained by tuning the relative size between the disks and the billiard. We define the volume 
V of the cavity and the pressure P, i.e., the derivative of each eigenvalue with respect to V. 
Reflecting the fact that the energy spectra of eigenvalues versus the disk size show a multitude of 
level repulsions, P — V characteristics shows the anomalous fluctuations accompanied by many van 
der Waals-like peaks in each of individual excited eigenstates taken as a quasi-equilibrium. For 
each eigenstate, we calculate the expectation values of the square distance between two disks, and 
point out their relationship with the pressure fluctuations. 
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I. INTRODUCTION 

The study of classical and quantum billiards has a long history [l|, I2J, |3|, y, |5( . From a 
viewpoint of classical dynamics, cavity billiards are classified into two categories according 
to their integrability or nonintegrability [6]. The stadium is a prototype of the nonintegrable 
and chaotic billiards, whereas highly-symmetric billiards like circle, triangle and square are 
integrable and regular. The quantum-mechanical feature of these billiards such as level 
statistics is nowadays well known [?], 

Most of the studies so far, however, are limited to the systems of a single point particle or a 



single disk confined in billiards. There exists 



number of mutually-interacting particles 







ittle work on the billiards which contain a finite 



Id . . The interaction between particles will 



make the system nonintegrable and chaotic even if the confining cavity is highly symmetric. 
On the other hand, recent development in hightechnology has fabricated quantum dots and 
optically-trapped atoms where a finite number of interacting particles are trapped in a small 
area. Quantum mechanics of these systems constitutes a topical subject 3, Q|- 

Considering the above circumstances, we want to understand what kind of novel phe- 
nomena should occur in the cavity billiards where interacting small particles are embedded. 
In this paper, we consider a pair of identical hard disks of finite size accommodated in a 
circular cavity. Each disk is assumed to bounce elastically with its partner and with the 
boundary wall, but otherwise showing a ballistic motion. Although the system is circularly 
symmetric, it is nonintegrable and chaotic thanks to the (short-range) interaction between 
the disks, as proved in the following analysis of classical dynamics. 

The tunable control parameter of the system is an aspect ratio of the radii between each 
disk and the circular cavity. When the aspect ratio is varied, how do the quantal eigenvalues 
and eigenfunctions behave? A multitude of level repulsions and level-spacing distributions 
like Wigner distribution will be shown to appear. Recalling the similarity between the 
present system and a nonideal gas confined within a container, the most important feature 
of this system would be a quantal equation of states, i.e., the volume V dependence of the 
pressure P at the cavity wall. Under a fixed radius of each disk, tuning of the aspect ratio 
corresponds to contraction or expansion of the volume V of the circular cavity. The pressure 
P at the wall boundary is obtained from the derivative of parameter-dependent eigenvalues 
with respect to V. In the classical nonideal gas theory, P does not decrease monotonically 
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as V is increased, and is accompanied by the van der Waals peak. This peak appears due 
to the competition between the finite size of each molecule and the attractive interaction 
between molecules. In the present system, we find the competition between the finite size 
of each disk and the quantum-mechanical correlation due to symmetrization of two-particle 
wavef unctions. Therefore P — V characteristics in the present system is expected to exhibit 
novel fluctuations unseen in ordinary quantum billiards that contain only a single particle 
or disk. This is an advantage of the interacting disk systems over the ordinary quantum 
billiards. By using two-particle wave functions, we shall evaluate expectation values for the 
distance between disks, which will elucidate a physical mechanism for fluctuations in the 
quantal equation of states. 

The organization of the paper is as follows: In Section |H1 we shall introduce a model 
system and analyze the underlying classical dynamics. We shall compute the maximum 
Lyapunov exponent with use of a modern technique by Gaspard and Beijeren [15.]. In Sec- 
tion we define the two-particle basis wavefunctions appropriate to quantize this compli- 
cated system and describe some technical details to obtain a standard form of the eigenvalue 
equation. The result for eigenvalues is given for various aspect ratios in Section IIVI Here 
we shall also analyze P — V characteristics, which will be the central theme of our study. 
Section |V| is devoted to investigation of expectation values for the distance between disks, 
where the single-particle density and two-particle correlation functions are exploited. The 
final Section is concerned with summary and discussions. 



II. MODEL SYSTEM AND CLASSICAL DYNAMICS 

Figure ^ illustrates our model system consisting of two identical hard disks moving within 
a circle billiard. Each disk shows a ballistic motion except when bouncing with its partner 
and with the billiard boundary. Let R and a(< R/2) be the radii for the circle billiard 
and each disk, respectively. Here the internal self-rotation (spinning motion) of each disk 
are ignored. The single-disk motion in the circle billiard is integrable since the number of 
degrees of freedom (two) accords with that of constants of motion, i.e., energy and angular 
momentum. However, the two-disk system proposed here has the degree of freedom (four) 
which is larger than the number of constants of motion, i.e., the total energy, angular 
momentum. Therefore the system becomes nonintegrable, and may be chaotic. 
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The chaotic dynamics is characterized by the positive maximum Lyapunov exponent 
A. This can be obtained by numerical simulation for a given initial condition by using 
linearized equations for deviations of an orbit. Although the degrees of freedom is four 
in the present system, the calculation of A is greatly simplified because of the nature of 
billiard motion [1J, [l5[. For the evaluation of A, each orbit is calculated over iV co i = 10 6 
collisions. Generally, A depends on the initial conditions. We randomly choose N- m n = 10 5 
initial conditions satisfying L z = 0. The average and standard deviation of obtained A's 
are shown in Fig. For < o = a/R < 0.4, the positive maximum Lyapunov exponent 
can be observed. The standard deviation systematically decreases below o ~ 0.25 as iV co i 
increases. The convergence, however, is not observed above a ~ 0.25 as shown by large 
standard deviations; This indicates that the ergodicity is not guaranteed in this region. 

Through the numerical simulation, we found that a small amount of initial conditions 
leads to A = indicating the existence of tori in the phase space in the full region of a. 
The ratio of tori to the whole phase space is, however, extremely small, and is estimated 
typically as smaller than of 10~ 3 . Although the present system should be said exactly to be 
a mixed system, the chaotic sea occupies almost whole part of the phase space, and one can 
proceed to investigate a quantum analogue of chaos in this system. 




FIG. 1: Model system. 
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FIG. 2: The ratio a = a/R vs. the maximum Lyapunov exponent A calculated by the numerical 
simulation. The total L z is taken as zero. The error bar means the standard deviation of obtained 
A's. 

III. METHODOLOGY OF QUANTIZATION 

To make variables dimensionless, we choose the scaling for coordinates and energy as 

r , a 



a 



R — a R — a' 

R „, E 

R — a 



x = E ' = — (!) 

2m(R-a) 2 

Note that by this scaling the ratio a = aj R = a'/R' is not changed. By the above scaling, 
the admissible range for r' = \r'\ becomes < r' < 1. Hamiltonian is then given in a 



dimensionless form as 



2^ { fi r >2 + r > f) r ' + r '2 fif)>2 j + ^ 

i=l,2 V » * * * * J i=l,2 



nir'r-r',1), (2) 



where 



and 




t/(r<) = ~ - " , (3) 



, (£ > 2a') 

V® = , ■ (4) 

oo (f < 2a!) 

The second and third terms in (J2J) represent confining potential by the hard wall of the 
circle billiard and inter-disk interaction due to their hard cores, respectively. In the case of 



a single disk inside the circle billiard, the appropriate wavefunction satisfying the boundary 
condition is given by Bessel function of integer order for the radial part multiplied by the 
angular function as 

J k {Xknr')e ike ', (5) 

where k and n are integers and Jfc(Afc n r') vanishes air' — 1, namely at the boundary wall. By 
choosing the wavefunction (jSJ) for each of the disks, Dirichlet boundary condition represented 
by the potential £7(r^) is automatically satisfied. On the other hand, the basis functions for 
the two-disk system are given by a product of (0) and take 

$ Q = J fcl (A fcini r' 1 )J fe2 (A fc2n2 ^)e jfcl ^e tfe2 ^/(X), (6) 

which are as yet not orthonormal. The multiplicative factor f{X), which represents the 
excluded volume effect caused by the hard cores of the disks, is defined by 

f 1 - e-P* (X > 0) 
f(X) = . (7) 

{ o (X < 0) 

Here X = r'^ + r'^ — 2r[r' 2 cos<p — Zq with <p = O^ — O^ and / = 2a'. X stands for a void between 
two hard disks, is a relative angle, and Iq is the admissible minimum distance between 
the centers of two hard disks, at which they touch each other. (3 is arbitrary positive real 
number, and may be called as the fictitious inverse temperature. We will discuss about (3 at 
the end of this Section. $ a vanishes at X = 0, satisfying another boundary condition that 
the wavefunction for two disks should vanish when they touch each other. In other words, 
the inclusion of the factor f(X) is approximately identical to incorporating the effect of the 
short-range repulsive interaction 1^ ( | — | ) in ©• Finally we symmetrize the wavefunction 
noting the indistinguishability of two disks. We here choose a symmetrized wave function 
by assuming the invariance of wavefunctions against the exchange of the disk coordinates. 
This choice is appropriate when the disks are Boson or Fermion forming a spin-singlet state. 
Therefore we use the basis function defined by 



^ a = (<l> a (l,2) + <l> a (2,l))/v / 2. (8) 
With use of (jHI) we shall proceed to construct the energy matrices. First, by operating 
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the Hamiltonian (J2J) on (jHJ), we find 

1 

+ 7f 



-8/3J fcl J fe2 e ifcl ^e ifc2 ^ e -^(1 - ^(r^ 2 + r^ 2 - 2/^ cos0)) 
_40 c iWi c <Wi c -*x (j^^-ft - r' 2 cos0) + Jj^te - rj cos, 



-2 J> . (9) 

Second, multiplying by "6ra" that corresponds to the complex conjugate of (jHJ), we 
integrate each term over coordinates. The integration about angles 9[, 6' 2 can be performed 
analytically. But we have to carry out the numerical integration about the radial coordinates 
r[, r' 2 . Third, noting that the total angular momentum L z is a good quantum number, we 
diagonalize the energy matrices for each of the fixed value L z . That is k% + k 2 = k[ + k' 2 = L z 
with L z = 0, ±1, ±2, and so on. 

Finally, before proceeding to diagonalization of each energy matrix, it must be regularized 
because we are using non-orthogonal basis functions. This is performed as follows: In 
the eigenvalue problem = E^> under consideration, we substitute the expansion ^ = 
C a ^a with the non-orthogonal basis functions. Then one reaches 

Hg'aCg — EN a t a C a , (10) 
a a 

With Ha! a, = (^a'\H\^ a ) and (^ a '\^a) = N a , a ^ Sofa. 

With use of the diagonalized norm-kernel Nj obtained by 

U~ l NU -> N d , (11) 
we introduce the regularized energy matrix 

H = N^U^HUN^ 2 . (12) 
As a consequence, we arrive at a standard form of the eigenvalue equation 

Hg = Eg, (13) 
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where g = N^U^C. This equation will be solved in the next Section. 

Before closing this Section, we should comment on the fictitious inverse temperature (3 
involved in the excluded- volume effect factor f(X) in (J7J). We find the calculated eigenvalues 
are stable against variation of (3 so long as it falls in 0.5 < (3 < 3.0 so that we shall fix to 
P = 1.5. 



IV. ENERGY SPECTRA AND PRESSURE 
A. Eigenvalues and level-spacing distribution 

We prescribe a = a/R with < a < 0.5 as a tunable parameter controlling the size of the 
disk. Under a fixed angular momentum L z , we diagonalize the energy matrix in 10 3 x 10 3 
dimensions by varying a . The lower half of 10 3 eigenvalues, which have a precision of 4 
digits, is used for our study below. In the manifold L z = 1, the energy spectrum against 
a is given in Fig. 3. The upward shifts of energies as a whole with increasing a reflect 
the confining of disks in a rapidly-decreasing effective area inside the wall boundary. This 
spectrum together with its partial magnification shows a multitude of avoided level crossings 
(level repulsions). Except for the lower energy region, the avoided crossings are widely seen 
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(a) (b) 




s 



FIG. 4: Histogram of level-space distribution: (a)Noninteracting disks with a = 0.2 and L z = 1; 
(b)Interacting disks with a = +0 and L z = 1; (c)Interacting disks with a = 0.2 and L z = 1; 
(d) Interacting disks with a = 0.2 and L z = 0. The solid, dashed and dotted curves show the 
Poisson distribution, Wigner distribution and the mixture of two independent Wigner distributions, 
respectively. 



in the full range of < a < 0.5. 

We proceed to investi gate the level-spacing distributions, where the usual unfolding pro- 
cedure is performed [bsl ll?)]. Figures 4 (a)-(c) are the level-spacing distributions P(s) for 
typical values of a with L z = 1, obtained from intermediate 450 levels by suppressing both 
the lowest 50 levels and the upper half of 10 3 levels. Figure 4 (a) is the result for the nonin- 
teracting two-disk system within the circle billiard, which approximately provides Poisson 
distribution P(s) = exp(— s). Small deviation from the Poisson distribution is however ob- 
served; We infer that this deviation comes from the particular properties of the zero points 
of the Bessel function which tends to be periodic in the asymptotic limit. In the case of 
interacting disks under consideration, the Wigner distribution P(s) = |exp(— its 2 ) can be 
seen in the full range of < o < 0.5. This distribution is known as a quantal signature 
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of chaos. It is interesting that Wigner distribution is obtained even in the point-disk limit, 
a = +0. This is due to the fact that, in this limit, the present system has a short-range 
interaction of delta-function type, and is distinct from a noninteracting two-disk system. 
The stability of level statistics against the variation of a is very convenient in our study 
below on the pressure. 

Before proceeding further, however, we should note an atypical level statistics in the 
manifold with L z = 0. Figure 4 (d), with a finite weight at the origin, is neither Wigner 
nor Poisson distribution. This distribution can be fitted by the mixture of two independent 
Wigner distributions (shown by the dotted curve in Fig. 4 (d)). We expect that the origin 
of this feature comes from the existence of a kind of the time-reversal symmetry. We can 
define a operator T commuting with the Hamiltonian as 

k,ni,— k,ri2 ^ k,ii2,— k,n\i 

where (k,n) denotes the one-particle state, and ^k, ni ,-k,n 2 is a (symmetrized) two-particle 
basis function whose total L z equals to zero. The eigenstates can be divided into two 
kinds by the eigenvalue of T, and the level spacing distribution becomes the mixture of two 
independent Wigner distributions. 

B. Pressure versus volume 

When a disk hit the cavity wall elastically, its momentum is reversed, which is compen- 
sated by the impulse on the wall. This is the origin of pressure on the wall. In the classical 
ideal gas theory, the pressure P is inversely proportional to the volume V. In the nonideal 
gas, on the other hand, there exists a competition between two kind of interactions, namely, 
the short-range repulsion due to hard cores of constituent molecules and the long-range at- 
tractive interaction between them. This competition leads to a peak called van der Waals 
peak in P — V characteristics. In the quantum two-disk system under consideration, we 
have both the hard-core repulsion and the strong quantum correlation due to symmetriza- 
tion of two-disk wavefunctions, and naturally can expect van der Waals-like peaks in P — V 
characteristics. The pressure P in each excited state taken as an equilibrium is obtained by 
means of the a-dependent eigenvalues. We choose a — 1, and define the cavity volume as 
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FIG. 5: Pressure versus cavity volume. From the bottom to the top, the level number is taken as 
1, 19, 38 and 45. 



V = ttR 2 . Then, the pressure is calculated as 

dE t dE t da dR 1 dE\ 
= ~~dV = ~~da~dRdV = 2ttR 3 da ' ^ 

where Ei stands for the energy for level I. Figure 5 shows P — V characteristics for several 
eigenstates. From the definition ()14)) . the van der Waals-like peaks and pressure fluctuations 
are attributed to level repulsions. In the low-lying states, there is no van der Waals-like peak, 
reflecting the absence of level repulsions. In the high-lying states, however, the number of 
peaks are increased, in contrast to the classical nonideal gas theory that accommodates only 
a single peak. 

A comment should be made here: In the case of a single disk in chaotic billiards like a 
stadium, one might also see the avoided crossing (AC) and level fluctuations as the system's 
parameter is varied. However, it is difficult to find its analogy to the nonideal gas because of 
the absence of inter-disk interactions. Furthermore, most of chaotic billiards have distorted 
wall boundary, for which there exists no global pressure as defined in (I14|) in the case of 
circle billiard. Thus P — V characteristics and the analogy to the nonideal gas should be 
meaningful only in the case of interacting particles (disks) confined in the circle billiard. 
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FIG. 6: (a) The pressure versus volume for a level 17 (the solid curve) and 18 (the dashed curve), 
(b) Corresponding average of the distance between two disks. The arrow shows the level crossing 
between these two levels. 

V. EIGENFUNCTION 



The fluctuations of the pressure mentioned above is a promising manifestation of quantum 
chaos in interacting two-disk systems in the billiard. In this section, we discuss this pressure 
fluctuation in terms of the wave function. 

Although the two disks originally have repulsive correlation due to hard-core interaction, 
quantum effective exchange interaction gives an additional strong quantum correlation be- 
tween disks. The averaged distance between disks is determined by the competition between 
repulsive hard-core interaction and quantum correlation depending on the control parame- 
ter a. We now wish to make clear whether this competition is responsible for the van der 
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Waals-like peak in quantum systems. For this purpose, we investigate the expectation values 
of the square distance 



for each level, where d — \r[ — r' 2 \. 

We concentrate on typical two neighboring levels employed in Fig. 6 (a), which bear sharp 
avoided crossings (AC) at the point indicated by the arrow. The avoided level crossing leads 
to the van-der- Waals-like peak (or dip). In order to study of the origin of this peak, we show 
the average of the distance between the disks in Fig 6 (b). Clearly, the sharp changes in 
the pressure curve corresponds to the change in the averaged distance; When the averaged 
distance increase (decrease) sharply, the pressure increases(decreases) similarly. This result 
is intuitively reasonable; If the averaged distance increases, the repulsive correlation between 
disks works strongly, and the existence probability of the disks near the circular wall will 
increase. Hence, the repulsive correlation contributes to the upward jump in the pressure 
curve. 

VI. SUMMARY AND DISCUSSIONS 

Quantum mechanics of a pair of identical hard disks confined in a circular billiard is 
investigated. Although the system is circularly symmetric, the short-range interaction be- 
tween the disks turns out making it nonintegrable and almost chaotic. With use of the 
two-particle basis wavefunctions that incorporate the excluded volume effect, eigenvalues 
and eigenfunctions are obtained for various ratios of the radii between each disk and the 
circle billiard. The energy spectra of eigenvalues versus the disk radius show a multitude of 
level repulsions. We find the Wigner level spacing distribution in the full range of the disk 
radius < a < R/2. The Wigner distribution in the vicinity of a = +0 should be ascribed 
to the strong inter-disk repulsion prevailing even at a = +0. 

The most important feature of this system is a quantal equation of states, i.e., the volume 
V dependence of the pressure P. The tuning of the aspect ratio corresponds to contraction 
or expansion of the cavity volume V. This tuning also changes the pressure P felt by 
the wall boundary. P — V characteristics shows novel fluctuations coming from quantum 
level repulsions. The van der Waals-like peaks in this fluctuations are attributed to the 




13 



competition between the finite size of two disks (responsible for the short-range repulsion) 
and the quantum-mechanical attractive correlation between them. The expectation values 
of the square distance between the disks is found to elucidate a mechanism for the pressure 
fluctuations. 
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